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Abstract 

FRW models of the universe have been studied in the cosmological 
theory based on Lyra's manifold. A new class of exact solutions has been 
obtained by considering a time dependent displacement field for variable 
deceleration parameter from which three models of the universe are de- 
rived (i) exponential (ii) polynomial and (iii) sinusoidal form respectively. 
The behaviour of these models of the universe are also discussed. Finally 
some possibilities of further problems and their investigations have been 
pointed out. 
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1 Introduction 

Einstein proposed his general theory of relativity, in which gravitation is de- 
scribed in terms of geometry, and it motivated the geometrization of other 
physical fields. One of the first attempts in this direction was made by Weyl [Tj 
who proposed a more general theory in which gravitation and electromagnetism 
is also described geometrically. However, this theory was never considered se- 
riously as it was based on the non-integrability of length transfer. Later Lyra 
[2] suggested a modification of Ricmannian geometry by introducing a gauge 
function which removes the non-integrability condition of the length of a vector 
under parallel transport. Subsequently, Sen [5] & Sen and Dunn [3] proposed 
a new scalar tensor theory of gravitation. They constructed an analog of the 
Einstein field equations based on Lyra's geometry which in normal gauge may 
be written as 

13 3 

R ij - 2 gi i R + ~ l9ij^>k<f> = -8nGT i: j, (1) 

where 4>i is the displacement vector and other symbols have their usual meaning 
as in Ricmannian geometry. 
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Halford [5] pointed out that the constant displacement vector field fa in 
Lyra's geometry plays the role of a cosmological constant in the normal general 
relativistic treatment. Halford [5] showed that the scalar-tensor treatment based 
on Lyra's geometry predicts the same effects, within observational limits, as in 
Einstein's theory. Several authors [7] have studied cosmological models based 
on Lyra's geometry with a constant displacement field vector. However, this 
restriction of the displacement field to be a constant is a coincidence and there 
is no a priorireason for it. Singh et al. [5] have studied Bianchi type I, III, 
Kantowski-Sachs and a new class of models with a time dependent displacement 
field and have made a comparative study of Robertson- Walker models with a 
constant deceleration parameter in Einstein's theory with a cosmological terms 
and in the cosmological theory based on Lyra's geometry. Recently Pradhan 
et al. [9] and Rahaman et al. [10] have studied cosmological models based on 
Lyra's geometry with constant and time time dependent displacement field in 
different context. 

The Einstein's field equations are a coupled system of highly nonlinear dif- 
ferential equations and we seek physical solutions to the field equations for their 
applications in cosmology and astrophysics. In order to solve the field equations 
we normally assume a form for the matter content or that space-time admits 
killing vector symmetries [TT] . Solutions to the field equations may also be gen- 
erated by applying a law of variation for Hubble's parameter which was proposed 
by Berman [T2| . In simplest case the Hubble law yields a constant value for the 
deceleration parameter. It is worth observing that most of the well-known mod- 
els of Einstein's theory and Brans-Dicke theory with curvature parameter k = 0, 
including inflationary models, are models with constant deceleration parameter. 
In earlier literature cosmological models with a constant deceleration parameter 
have been studied by several authors [9] Q3] ■ But redshift magnitude test has 
had a chequered history. During the 1960s and the 1970s, it was used to draw 
very categorical conclusions. The deceleration parameter go was then claimed 
to lie between and 1 and thus it was claimed that the universe is decelerating. 
Today's situation, we feel, is hardly different. Observations (Knop et al. [Tl] ; 
Riess et al. [TS]) of Type la Supernovae (SNe) allow to probe the expansion 
history of the universe. The main conclusion of these observations is that the 
expansion of the universe is accelerating. So we can consider the cosmologi- 
cal models with variable cosmological term and deceleration parameter. The 
readers are advised to see the papers by Vishwakarma and Narlikar [16] and 
Virey et al. [17j and references therein for a review on the determination of the 
deceleration parameter from Supernovae data. 

Recently Pradhan et al. [18] have studied the universe with time dependent 
deceleration parameter in presence of perfect fluid. Motivated by the recent re- 
sults on the BOOMERANG experiment on Cosmic Microwave Background Ra- 
diation (Bernardis [10]), we wish to study a spatially flat cosmological model. In 
this paper, we have investigated spatially non-flat and flat cosmological mod- 
els with a time dependent displacement field within the framework of Lyra's 
geometry. We have obtained exact solutions of the field equations of Sen [3] 
by taking the deceleration parameter to be variable. This paper is organized 
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as follows. The metric and the field equations are presented in Section 2. In 
Section 3 we deal with a general solution. The Sections 4, 5, and 6 deal with the 
three different cases for the solutions in exponential, polynomial and sinusoidal 
forms respectively. Finally in Section 7 concluding remarks arc given. 



2 Field Equations 

The time-like displacement vector fa in the equation ([1]) is given by 

fa = (0,0,0, f3(t))- (2) 
The energy-momentum tensor in the presence of a perfect fluid has the form 

Tij = (p + p)utUj - pg i:j . (3) 

1, where u f = (0,0,0, 1). The 



together with co-moving coordinates u % Ui 
metric for FRW spacetime is 



ds 2 = dt 2 - R 2 (t) 



dr 2 



1-kr 2 



\de 2 



sur 



(4) 



where, k = 1, —1,0. For this metric, the field equations (JTJ) with the equations 
HJ and flSJ) take the form 



, tt2 3k 3(3 2 

3# + T^T - — 



R 2 



2H + 3H 2 



k 

^2 



3^ 
4 



-XP, 



(5) 
(6) 



where \ = 8ttG and H = R/R is the Hubble's function. Equations (f5j) and 
lead to the continuity equation 



X p+ -13/3 + 3 



3 

X(P + P) + 



H = 0. 



Assuming an equation of state 

p = 7/5, < 7 < 1 
Eliminating p(t) from and © we obtain 

2H + 3(1 + 7 )ff 2 + (1 + 3 7 ) A + |(i _ 7)/3 2 



0. 



(7) 



(8) 



(9) 



Here (3 2 plays the role of a variable cosmological term A(£). We have two 
independent equations in three unknown viz R(t), p(t) and (3. Therefore we 
need one more relation among the variables in order to obtain a unique solution. 
Hence, we consider the deceleration parameter to be time dependent. 
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3 Solutions of the field equations 

We consider the deceleration parameter to be variable 



RR I H + H 2 



Q = -— = -[ J T2 )=b (variable). (10) 



The equation (|10p may be rewritten as 

R R 2 

5 + ^-°- ("J 

The general solution of Eq. (fTTj) is given by 

J i dR dR = t + to, (12) 

where m is an integrating constant. 

In order to solve the problem completely, we have to choose / ^dR in such 
a manner that Eq. (fT2")) be integrable. 

Without any loss of generality we consider 

J ^dR = In L(R), (13) 

which does not effect the nature of generality of solution. Hence from Eqs. Q12[) 
and (fP3")) . we obtain 

J L(R)dR = t + m. (14) 

Of course the choice of L(R), in Eq. (|14[) . is quite arbitrary but, since we are 
looking for physically viable models of the universe consistent with observations, 
we consider the following case: 

4 Solution in the Exponential Form 

Let us consider L(R) = where k\ is arbitrary constant. 

In this case on integration of Eq. (fH)) gives the exact solution 

R{t) = k 2 e kl \ (15) 

where k 2 is an arbitrary constant. Using Eqs. (JSJ) and (fT4|) in ([9]) and §5§ or ([6]), 
we obtain expressions for displacement field /3, pressure p and energy density p 
as 

2 _ 4(l + 7 )fc^ 4(l + 3 7 )fc , . 

P (l- 7 ) 3(l- 7 )fc|e 2fc i t ' 1 j 



27 

XP = X7P 



(1-7) 



3k\ 



k 2 e 2k ^ 



(17) 
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Figure 1: The plot of scale factor R(t) vs time with parameters k\ = 0.01, 
k 2 = 0.5, and 7 = 0.5 

From Eq. (fT5"|) . since scale factor can not be negative, we find R(t) is positive 
if k 2 > 0. From Figure 1, it can be seen that in the early stages of the universe, 
i.e., near t = 0, the scale factor of the universe had been approximately con- 
stant and had increased very slowly. At specific time the universe had exploded 
suddenly and expanded to large scale. This is consistent with Big Bang scenario. 

From Eq. (fTl)]) , it is observed that /3 2 is a decreasing function of time. As 
mentioned earlier the constant vector displacement field 4>i in Lyra's geome- 
try plays the role of cosmological constant A in the normal general relativistic 
treatment and the scalar-tensor treatment based on Lyra's geometry predicts 
the same effects, within observational limits, as the Einstein's theory. Recent 
cosmological observations (Garnavich et al. [50] , Perlmutter et al. [2T], Riess 
et al. [22], Schmidt et al. [23]) suggest the existence of a positive cosmologi- 
cal constant A with the magnitude A(Gh/c 3 w 10 -123 . These observations on 
magnitude and red-shift of type la supernova suggest that our universe may be 
an accelerating one with induced cosmological density through the cosmological 
A-term. In our model, it is seen that (3 plays the same role as cosmological 
constant and preserves the same character as A-term, in turn with recent ob- 
servations. 

From Eq. (fT7]) . we observe that p > and p > for k > 0. We also 
see that the energy density decreases to a small positive value and remains 
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constant thereafter. The expressions for 1 and p cannot be determined for the 
stiff matter (p = p) models. The expressions for f3 2 and p corresponding to 
7 = 0, 1/3 are given in Table 1. 
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Table 1: Values of /3 2 and p for dust and radiation exponential models. 

We can obtain the values of 1 and p for flat FRW model if we set k = 
in Eqs. (JTHJ) and fTTfl . 



5 Solution in the Polynomial Form 



Let = — j= i , , where fcq and fei are constants. 

In this case, on integrating, Eq. ([H]) gives the exact solution 

i?(t) = ait 2 + a 2 t + Q'3, 



(18) 



where a.\, a 2 and 0:3 are arbitrary constants. Using Eqs. ((5J) and l|18p in © 
and ([5]) or we obtain the expressions for displacement field f3, pressure p 
and energy density p as 



ltoijag + (2 + 37)ait 2 + (2 + 3 7 )q 2 £] + 4(1 + 3t)(q| + fc) 
3(7- l)(ait 2 + a 2 i + a3) 2 

4 7 [5a 2 t 2 + 5aia 2 < + a| + ai«3 + fe] 



XP = X7P 



(l-7)(ait 2 + a 2 i + Q!3) 2 



(19) 
(20) 



From Eq. (fT8)) . we note that R(t) > for < t < 00 if a\, a% and «3 are 
positive constants. Figure 2, shows that the scale factor is a decreasing function 
of time, implying that our universe is expanding. 



Eq. (|19|) shows that (3 2 < for all times as 7 — 1 < and is a decreasing func- 
tion of time, characteristically similar to A in Einstein's theory of gravitation. In 
this model, (3 plays the role as cosmological constant and it preserves the same 
character as A-term. This is consistent with recent observations (Garnavich ct 
al. [20], Perlmutter et al. [21], Riess et al. [22], Schmidt et al. [23]). A negative 
cosmological constant adds to the attractive gravity of matter; therefore, uni- 
verse with a negative cosmological constant is invariably doomed to re-collapse. 
A positive cosmological constant resists the attractive gravity of matter due to 
its negative pressure. For most of the time, the positive cosmological constant 
eventually dominates over the attraction of matter and drives the universe to 
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Figure 2: The plot of scale factor R(t) vs time with parameters a\ — 1.00, 
a 2 = 4.00, a 3 = 1.00 and 7 = 0.5 



expand exponentially. 

The expressions for 1 and p cannot be determined for the stiff matter 
(p = p) models. For dust model (7 = 0), 1 and p(t) are given by 

2 lgoi [2ait 2 + 2a 2 t + q 3 ] + 4o| + 4/c 
^ " (ait 2 + a 2 * + a 3 ) 2 ' 1 j 



4[5a 2 i 2 + 5aia 2 t + a\ + aia$ + k 
(ait 2 + a 2 t + a 3 ) 

For radiative model (7 = g), /3 2 and p(t) are given by 



2 _ 801 [2 ai t 2 + 2a 2 t + q 3 ] + 8o| + 8fc 
; ~ (a^ 2 +a 2 t + a 3 ) 2 ' 1 ' 



6[5a 2 t 2 + ha\a 2 t + a 2 + a\a^ + k] 
(ait 2 + a 2 t + a 3 y 

If we set k = 0, in above equations (TIT)!) " P^l - we get solutions for flat FRW 



™ = ~ ^X^7T^2 ( 24 ) 



universe. 
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6 Solution in the Sinusoidal Form 

Let L(R) = -—=5, where (3 is constant. 

In this case, on integrating, Eq. (fT4|) gives the exact solution 

R = Msm(/3t) +Ncos(pt) + j3 u 



(25) 



where M, N and f3± are constants. Using Eqs. © and (|25p in © and ((SJ) or 
((HJ) , we obtain the expressions for displacement field f3, pressure p and energy 
density p as 



2 



where 



A[2(3 2 (M 2 + N 2 ) + 20 2 1 (P - ft) - 3(1 + 7 )Q - (1 + 3 7 )fc] 
3(1 - 7 )P 2 

2 7 [3Q + 2fc - /3 2 (M 2 + N 2 ) - f3 2 1 {P - fa)] 



XV = XI P 



(1- 7 )P2 



(26) 



(27) 



P = M sin ftt + N cos fo + fa, 
Q = (M cos (3t - N sin (3t) 2 . 
From the Figure 3, we note that at early stage of the universe, the scale of 




Figure 3: The plot of scale factor R(t) vs time with parameters M = 2.00, 
N = 1.00, (3 = 10.00, /3i = 0.2, and 7 = 0.5 

the universe increases gently and then decreases sharply, and after wards it will 
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oscillate for ever. We must mention here that the oscillation takes place in pos- 
itive quadrant. This has physical meaning. 



The expressions for 1 and p cannot be determined for the stiff matter 
(p = p) models. For dust model (7 = 0), (3 2 and p(t) are given by 

R 2 _ 4[2/3 2 (Af 2 + N 2 ) + 2(3 2 (3 1 {P - (3,) - 3Q - k] 

P ^2 . (28) 

2[3Q + 2k - (3 2 (M 2 + N 2 ) - 0*fk(P - ft)] 
XP= pa ' ( 29 ) 

For radiative model (7 = |), /3 2 and p(t) are given by 

p2 = 4[/3 2 (A/ 2 + N 2 ) + [Pfh(P - A) ~ 2Q - fc] (30) 

_ 3[3Q + 2k- (3 2 (M 2 + TV 2 ) - P 2 P,{P - A)] 

XP p 2 1 \ 6L ) 

For flat FRW universe, we put k = in above results. Since, in these cases, we 
have many alternatives for choosing values of M, N, (3, /3i, it is just enough to 
look for suitable values of these parameters, such that the physical initial and 
boundary conditions are satisfied. 



7 Conclusions 

In this paper we have obtained exact solutions of Sen's equations in Lyra geom- 
etry for time dependent deceleration parameter in FRW spacetime. The nature 
of the displacement field (5{t) and the energy density p(t) have been examined 
for three cases (i) exponential form (ii) polynomial form and (iii) sinusoidal 
form. The solutions obtained in Sections 4, 5 and 6 are to our knowledge quit 
new. Here it is found that the displacement field plays the role of a variable 
cosmological term A. 

Recently there is an upsurge of interest in scalar fields in general relativity 
and alternative theories of gravitation in the context of inflationary cosmology 
(La and Steinhardt [24] ; Ellis [25]; Barrow [26]). Therefore the study of cos- 
mological models in Lyra geometry may be relevant for inflationary models. 
Further the space dependence of the displacement field (3 is important for inho- 
mogencous models for the early stages of the evolution of the universe. Besides, 
the implication of Lyra's geometry for astrophysical interesting bodies is still an 
open questions. The problem of equation of motion and gravitational radiation 
need investigation. Finally in spite of very good possibility of Lyra's geometry 
to provide a theoretical foundation for relativistic gravitation, astrophysics and 
cosmology, the experimental point is yet to be undertaken. But still the theory 
needs a fair trial. 
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